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Abstract. Redshift-space distortions are generally considered in the plane parallel limit, 
where the angular separation between the two sources can be neglected. Given that galaxy 
catalogues now cover large fractions of the sky, it becomes necessary to consider them in a 
formalism which takes into account the wide angle separations. In this article we derive an 
operational formula for the matter correlators in the Newtonian limit to be used in actual 
data sets. In order to describe the geometrical nature of the wide angle RSD effect on 
Fourier space, we extend the formalism developed in configuration space to Fourier space 
without relying on a plane-parallel approximation, but under the extra assumption of no 
bias evolution. We then recover the plane-parallel limit not only in configuration space 
where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections 
that should be included in large surveys as a perturbative expansion over the plane-parallel 
results. We finally compare our results to existing literature, and show explicitly how they 
are related. 
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1 Introduction 

Redshift-space distortion (RSD) is a general denomination for the corrections on distances 
attributed to objects from their redshift measurements, due to their peculiar velocities. Mod¬ 
els for the peculiar velocities rely on the theory of structure formation, in linear or nonlinear 
regimes. We shall here be interested in the linear regime, and therefore we will assume 
that dark matter distribution is well described by linear fluid equations, and that galaxies 
overdensities are associated to dark matter overdensities through a linear bias. 

Under these hypothesis, redshift space distortions become a geometrical effect since the 
corrected density field 6^ is related to the underlying linear density field 6 thanks to the 
integro-differential operator composed of a non-local part (the inverse of the Laplacian) and 
the radial part of the Laplacian 

(5^(r) = (5(r) ^dr[r‘^dr{A~^6{r))] . (1.1) 

When we study correlation functions, we derive from the properties of the operator in eq. (1.1) 
that the RSD corrections to the two-point function only depend on the shape of the triangle 
formed by the observer and the two objects being correlated. The limit in which the triangle 
is squeezed and close to isosceles corresponds to the plane-parallel approximation, or flat-sky 
approximation. The plane-parallel approximation is well justihed if the distance between the 
pair of galaxies is much smaller than the distances from the observer to each of the galaxies. 

In the plane-parallel limit, eq. (1.1) can be simplified and leads to the classical Kaiser 
formulation of the problem [2, 3]. 

Gurrent and future galaxy surveys cover sizable fractions of the sky such as BOSS [4] 
and Euclid [5]. For these catalogues, one cannot assume anymore that the separation angle 
between two galaxies is small, what makes the plane-parallel approximation unjustified and 
compels the incorporation of the so-called wide angle effects. Based on the analysis of simu¬ 
lated samples, [6, 7] have shown that, although sizable, wide angle effects may be of the size 
of systematic errors in configuration survey analysis if suitable geometrical configurations 
and analysis procedure are adopted. Hence, for practical purposes, the Kaiser formula would 
still provide accurate prescription for the RSD corrections. One of our goals in this paper 
is to build a perturbative expansion around the plane-parallel limit and show that the first 
corrections due to wide angles are quadratic in the ratio of the distance between the galaxies 
to the mean distance from the observer. 

The study of wide angle effects on RSD is not recent. A non-exhaustive list of works 
dedicated to the problem includes [8-14]. Building on those earlier articles, [1] was able to 
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calculate in full generality two-point correlation function, allowing to estimate the difference 
with the standard plane-parallel limit. 

Using the operator in eq. (1.1), the usual treatment of wide angle effects, as in [15], 
is to decompose the distorted matter density field (i^(r) in spherical harmonics, e.g. as 
for the Cosmic Microwave Background analysis. For every distance, this leads to define 
the associated multipoles since the redshift-space distortions do not break the 

statistical isotropy around the observer, the statistics is encoded in the diagonal part of the 
multipoles correlations, that is 

= (5££'-5wC'|(ri, r2). (1.2) 

Both the expression found with the conhguration space approach in [1] and the one 
derived with a multipole approach in [15] are much more complicated than the simple original 
Kaiser formula that holds in the plane-parallel limit. In this article, our second goal is to 
relate both approaches and to link them with our perturbative expansion around the plane- 
parallel limit. 

2 Overview of results 

2.1 Wide angle effects in confignration space 

The structure of the correlation function depends only on the shape of a triangle formed 
by the observer and the sources at ri and r 2 . The plane-parallel limit corresponds to a 
squeezed configuration of this triangle where r = r 2 — ri is much smaller in norm than ri and 
r 2 . In order to expand the correlation functions around this plane-parallel approximation, 
it appears more appropriate to consider the correlation functions as depending on r and the 
median distance d as illustrated in hgure 1, with the dehnitions 

e^(d, r) ^ C^(ri, r 2 ) = (<5^(ri)<5^(r2)) , r = ra - ri , d = . (2.1) 

In this article we give the general expression for the correlation function in configuration 
space, checking that we recover the results of previously existing literature, and we then 
expand it around its plane-parallel limit so as to grasp the structure of the wide angle 
corrections. We also show that instead of using the median position d as an average position, 
it is possible to use the bisector to dehne another type of average position, as it leads to the 
same plane-parallel limit. 

If the correlation function was statistically homogeneous, then it would depend only 
on r, and not on d. And if it was also statistically isotropic, it would actually depend 
only on one degree of freedom, r. This would be the case if RSD effects were ignored. 
However the distorted field is not homogeneous as it also depends on the velocity of the source 
with respect to the observer, and not just on the velocity independently. Nevertheless, the 
global rotational invariance around the observer removes three degrees of freedom, implying 
that the correlation function is only a function of three degrees of freedom which are r, d 
and //dr = h • d. In the plane-parallel limit, it depends on r and //, but not on d. It is 
thus appropriate to expand the general two point correlation function as a general angular 
multipole expansion 

OO OO 

^""(dji*) = ^ , with //dr = r-d. (2.2) 

n=0 i=0 
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The are the lowest order coefficients which arise in the plane-parallel limit, and to be 
more precise, only are non-vanishing. The describe then corrections due 

to wide angle effects. The natural small parameter for this expansion is the ratio between the 
distance between two points being correlated, and their average distance from the observer, 
that is r/d. The further the two points are with respect to the observer, the smaller the 
corrections are. The geometrical structure at each order is described by the angular variable 
which depends only on the relative directions between the difference of positions r and the 
average direction d, that is it depends only on fj.. We find that the coefficients Q ^ in 
the expansion (2.2) are non-vanishing only if I and n are either both odd or both even. 
Furthermore, we show that for the median and the bisector parametrizations of the average 
distance, these coefficients do not vanish only if i and n are both even. As a consequence, 
the first order corrections vanish, and these choices should thus be preferred to minimize 
the wide angle effects. In these cases, we compute explicitly the Q up to second order 
providing the first set of corrections to the plane-parallel limit. These coefficients emerge 
from the analysis of eq (3.21), explicitely written in eq (4.7), which constitutes one of the 
main results of this paper. 



Figure 1. Effect of Fourier transform on variables constrained to a triangular geometry. The Fourier 
conjugate of ri is ki, the Fourier conjugate of r2 is k2; the Fourier conjugate of d is Ak = ki — k2, 
while the conjugate of r is k. 


2.2 Wide angle effects in Fourier space 

In order to understand how such expansion arises in Fourier space, one must first realize that 
when performing a double Fourier transformation on a function whose variables are related in 
a triangular conhguration, the Fourier conjugate to the difference of positions r is the median 
Fourier modes k, and the Fourier conjugate to the median distance d is the difference of the 
Fourier modes Ak. This geometry of the Fourier space is illustrated in figure 1. This basic 
geometrical relation can be understood in two steps. First, for a > 0, the Fourier transform 
of C'^(ari,r 2 ) is ^C'^(^,k 2 ), i.e., dilatations (resp. contractions) in conhguration space lead 
to contractions (resp. dilatations) in Fourier space, and therefore if ri is shorter than r 2 , ki 
will be longer than k 2 - Secondly, when performing a Fourier transformation on a two point 
correlation function, we introduce the product 

with k = (ki -|- k 2 ) /2 , 
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Ak = ki - ka , (2.3) 






from where the crossed conjugate relation among median and difference of modes can be 
inferred. It is thus natural to dehne the correlation function in Fourier space using these 
variables as 

C^(Ak,k) = C^(ki,k2) = (5^(ki)<i-(k2)). (2.4) 

The homogeneity in Fourier space is expressed by the fact that the correlation function de¬ 
pends only on the average Fourier mode, and not on the difference. The inhomogeneity 
introduced by the RSD effects translates into the fact that in Fourier space there are off- 
diagonal correlations. In this article, we derive the general expression for the correlation 
function in Fourier space and exhibit the off-diagonal contributions. Following the corre¬ 
spondence (2.3), the corrections introduced should be expressed as an expansion in |Ak|/A: 
in the form 

1 f I /\ H \ ^ 

C(Ak,k)=MAk)ci">(m---™j;(!^ (2.5) 

' ' n=0 ^ ^ £=0 

(£,n)^(0,0) 

with /XkA = k • Ak. The correspond to the homogeneous contribution of the plane- 

parallel limit, for which only C^\ are non vanishing. The are the wide angle 

corrections which break homogeneity. For each order, the geometrical dependence is only a 
function of the angle between the average Fourier modes and the difference of the Fourier 
modes, and this is understood from the correspondence (2.3). We do not perform such 
expansion explicitly except for the lowest order corresponding to the plane-parallel limit. 


2.3 Wide-angle effects in mixed configuration/Fourier space 

In fact the RSD effects can also be apprehended using a mixed space, where the median 
distance is looked at in conhguration space, but the dependence in the separation of the 
sources is considered in Fourier space. This can be obtained either by Fourier transforming 
the r dependence, that is considering ^^(d,k) instead of in ,^^(d,r), or by inverse Fourier 
transforming the Ak dependence, that is by considering (^^(d,k) instead of C^(Ak,k). We 
check that both approaches lead to the same result as they ought to. In this mixed space 
the natural expansion is 


r(d,k) = c(<i,k) = j; 

n=0 


1 

kd 


n oo 


VRf)(A))P,(/rkc 


1=0 


/ikd = k • d , 


( 2 . 6 ) 


The functions V^^\k) are related to the configuration space coefficients Ci^\r) through an 
integral on spherical Bessel functions as 


J <irr'^{kr)'^ji{kr)^^''\r), (2.7) 

implying that they do not vanish only if i and n are either both odd or both even. Further¬ 
more, they are directly related to the {k) by numerical factors. Defining s by ^ = 2-|-n-|-2s, 
these relations are simply 


Vi'"\k) 


(27r)3/2 


2”+*(n-hs)! 
(2s+ 3)!! 


Ch(k) 


if (£, n) 7 ^ (0,0), 




1 ^^(0) 
(27r)3/2‘’0 


{k). 

( 2 . 8 ) 
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In the plane-parallel limit n = 0, and only V^\k), V^\k), V^\k) are non-vanishing, 
as they correspond to the standard Kaiser formula which we show to be only meaningfully 
dehned in this mixed space. Using eq. (2.7) we are able to deduce the coefficients v]p\k) 
of the mixed space expansion from the configuration space correlation function expansion, 
that is from the ^^^\r). There is then no need to compute the coefficients c!f'\k) of the full 
Fourier space expansion since we find that they can be easily deduced from the mixed space 
expansion thanks to eq. (2.8). With this mixed space approach, one is led to consider the 
various power spectra at each given distance d so as to be able to take into account wide 
angle effects, and the encode the angular dependence of the corrections at each 

order. In practice, the inclusion of the wide angle effects in data analysis implies that we 
must build a spectrum for each subset of pairs having a common average distance d. 

2.4 Outline of the article 

The paper is structured as follows: we start by deriving eq. (1.1), and use it to obtain the 
general form of the two-point correlation function. We show how to obtain the plane-parallel 
limit of the general expression and proceed with the analysis of its counterpart in Fourier 
space. The mixed configuration/Fourier spaces emerge naturally from the structure of the 
problem and we show that the plane-parallel limit is well defined in the mixed space, from 
where the Kaiser formula is consistently derived. We also investigate the decomposition of 
the distorted density fields in spherical harmonics and the Cfs both in configuration and 
Fourier spaces. This constitute a first part of the paper, where we show different ways to 
approach the problem, and the consistency among them. In the second part of the paper we 
construct a perturbative treatment of the wide angle effects over the plane-parallel results. 
We address the problem first in configuration space using the two-point correlation function 
and we exhibit the wide angle corrections using different parametrizations for the triangular 
made by the observer and the pair of galaxies. These results are then translated into mixed 
configuration/Fourier space, where we introduce the idea of spectrum at a given distance. We 
advocate that it is the most meaningful way to extend the Kaiser formula so as to incorporate 
wide angle effects. We also show that the general correlators in full Fourier space can be well 
dehned. 

As we shall see, although the wide angle problem can be treated in full generality in 
conhguration space, the complexity introduced by bias evolution clouds the description in 
Fourier space. We shall, therefore, only study Fourier and mixed descriptions under the 
simplifying hypothesis of no bias evolution, what limits the practical applicability of our 
work. It allows, nevertheless, the investigation of the intricate geometrical nature of RSD in 
Fourier space, which is one of our goals in this paper. 

3 General formalism of redshift-space distortions 
3.1 Physical origin 

We consider a perturbed Friedmann-Lemaitre (FL) universe with only scalar perturbations 
<h which is identihed with the gravitational potential. The corresponding metric is of the 
form 

ds^ = —dt^(l -|- 24>) -|- a(t)^(l — 2^)6ijdx^dx^ , (3.1) 

where a{t) is the scale factor. We then define the Hubble expansion rate as H = d/a, and 
we also introduce the conformal time defined by adrj = dt. The conformal Hubble rate is 
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then simply given hy % = aH. In a simple Newtonian description we ignore the effect of 
the gravitational potential <I> on the propagation of light and the redshifting of energies, but 
we consider its effect on a source velocity and the corresponding associated redshift. If the 
velocity of an emitting source is decomposed as 


= (l,u*) 


V = 


dx* 

drj 


then the redshift at which this source appears today ( 770 ) is given by 


I 


a(%) 


= va 


_ -(1 + ViX") , X* = — X = VX,;X‘ 

1 + z 0(770) X 

where a{i]e) is the scale factor at emission. Given that 


a(r 7 e)(l + UjX*) ~ a 


Ve + 


ViX" 


-HiVe 


(3.2) 


(3.3) 


(3.4) 


ignoring the effect of the peculiar velocity amounts to an error in the time of emission. This 
translates into an error on the conformal distance estimated from the background geometry 
through |xo — Xe| = t/o — Noting x = (x^,x^,x^) the true coordinates of an object at 
emission, and s = (s^, s^, s^) its inferred position when the peculiar velocity is ignored, these 
positions are related at hrst order in the velocity by 


s* = X* + a*, a* = q;,.x* , ar = —{vj.x^) = ^ . (3-5) 

Txe Txe 

If we were to include all the gravitational effects, then this relation would involve more terms, 
and a detailed analysis can be found in e.g. [15-19]. 

Since we consider only scalar perturbations, the velocity divergence 9 is characterized 
fully by the velocity field V through 


u* = d% 


div^ _ AG 


(3.6) 


The relation between V and the matter density is obtained from the resolution of Euler, 
continuity and the Poisson equations [20]. We then hnd that the velocity divergence 9 is 
directly related to the comoving matter density contrast S = p/p— 1, where p and p are 
respectively the background and the full energy density of matter. Indeed, for each Fourier 
mode we can relate them by a relation of the form 


9{k,t) = —f3{t)6{k,t) V{k,t) = . (3-7) 

The matter density p^ inferred from the redshifted coordinates s, is different from the underly¬ 
ing matter density p, and we thus want to relate the associated density contrast 5^ = {p^—p/)p 
to the true density contrast <5. The number of galaxies cannot depend on our coordinates 
system and the two matter densities are thus related by 


/7^(s)d^s = /3(x)d^x 


(3.8) 


Using the Jacobian of the transformation (3.5) which is d^s ~ 


(1 -|- (9ia*)d^x, we obtain 


p^{s) ~ p[x(s)][l - diO^] ~ p{s) - pdia\s) 


5^ = 5- diot 


(3.9) 
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since /o[x(s)] — p{s) is a second order quantity as the background matter density is homo¬ 
geneous. This means that at lowest order, the volume are affected by the divergence of the 
displacement a*, and this translates into a modification of the density with an opposite sign. 

Since the displacement field a* is only in the radial direction, it is useful to recast its 
divergence in spherical coordinates to get 

diO^ = ^drir'^ar) = :r^dr{r‘^drV) = --^dr[r^dr{A~^6)] (3.10) 

Her 

where we use the usual notation r = |r|, and therefore 

5^(r) = (5(r) ^dr[r‘^dr{A~^5{v))] . (3.11) 

In order to express the radial derivatives in this expression, we decompose the Laplacian in 
three dimensions into a radial operator and the Laplacian on the sphere according to 


A = Or + 


L 




= -;^drr 


^dr 


L? = l2 + 


_ f) 

1 _ 1/2 > 


lI = d^{l - 


(3.12) 

where u is the cosine of the polar angle and (p the azimuthal angle of f in a spherical 
coordinates system. When it brings no ambiguity, we will omit the index in L?. We then get 
from eq. (3.11) an equivalent relation between the two types of densities which is 


5^{r) = (1 + ^),5(r) - A-^6{r). (3.13) 

The time dependence of /3 can be written as a dependence on r thanks to the background 
geodesic relation r(t) which can be inverted into a relation t{r). Hence, in eq. (3.13) (3 
stands for j3{r) = I3{t{r)). The corrections considered in eq. (3.13) constitute a subset of 
all the relativistic effects in first order of perturbation theory. As observed by [21] - based 
on the analysis presented in [15] -, only density fluctuation, redshift space distortion and 
lensing can contribute significantly to corrections of number counts in the context that we 
are considering. Lensing, however, has a distinct geometrical nature and we shall not be 
discussed here. For simplicity we also assume a constant selection function. 


3.2 Redshift-space distortions in confignration space 
3.2.1 General expression of the correlation function 

We first start by deriving the two-point correlation function of distorted densities in configu¬ 
ration space, using the previous relation (3.11). We do not assume that the angle separation 
is small to begin with, so as to obtain the most general expression. In the next section, we 
then detail how we can recover the usual plane-parallel limit, known as the Kaiser formula. 
The statistics of the fundamental field 5 is known in Fourier space, as it is given by the matter 
power spectrum defined by 

C'(ki,k2) = (<5(ki)5*(k2)) = V{k)6D{\^i - k2), (3.14) 

where the 5(k) are the Fourier transform of 6 today. They are related to the Fourier modes 
at an earlier time by a transfer function as 

(5(k,t) = G(t)(I(k), G{to) = l, (3.15) 
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since when radiation can be neglected the modes evolve at the same pace [21]. We first 
need to obtain the distorted field 5^ in configuration space as a function of the underlying 
density field 5 in Fourier space. Given that the Fourier transform of is —(5(k)/A:^, we 

get immediately from eq. (3.11) that this relation takes the form 


<5^(r) = G(r) 


d^k 

(27r)^/^ 


m 




e 


ikr 


(3.16) 


where G{r) = G{t{r)). Correlating the distorted field in two different points, we use eq. (3.14) 
to remove one of the two Fourier integrals, and the other one is performed easily in spheri¬ 
cal coordinates once the exponential is expanded in spherical harmonics with the Rayleigh 
expansion (A. 8). We finally obtain 


C^(ri,r2) 



1-^0 

A;2 


joikr ), 


(3.17) 


where we used the notation j3i = /3(ri), (32 = /3(?'2) and Gi = G{ri), G 2 = G{r 2 ), and we 
have defined the difference of positions and the associated norm as 


r = r 2 -ri, = rf + - 2rir2Ui2 , 1^12 = ri • h . (3.18) 


In order to express the radial operators and Or 2 , we hrst use that 


Li. 


L? 


A = a. + -^ = Or, + . 


(3.19) 


We then note that the radial operators are only applied on a function of r in eq. (3.17), where 
r, ri and r 2 are in a triangular configuration. For any function /(r) 


OrJ{r) = A/(r) - ^L? /(r) = a/(r) - 


^/( 0 > Or2fir)=Orfir)-^^f{r) 
rf r| 


(3.20) 

where we have used i) that for a function of r only, the Laplacian reduces to Or, the partial 
derivative dr becoming in that case total derivatives, and ii) that the angular operators 
and do not depend on the azimuthal angle, as they depend only on the polar angle, and 
thus both reduce to • 

With these relations and eq. (A.7) we can recast (3.17) as 


C'^(ri,r2) = GiG2 


I IMkr)] = 


l + /3i + 


k'^dk 

2tt‘^ 

/3iL2 


V{k)I[jo{kr)] 


(kriY^ 


1 + /32 + 




{kr2Y _ 


jo{kr) 


(3.21) 


Note that this result could also have been obtained by using the relation (3.13) between the 
distorted and the underlying matter densities. This expression for the correlation function in 
configuration space is the first major result of this article. Indeed, it is from this expression 
that we will extract the plane-parallel approximation and its corrections due to large angle 
effects. 
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3.2.2 Plane-parallel limit 

In order to have a better insight on the general expression (3.21), we check in this section 
that we can obtain from it the standard result in the plane-parallel limit. 

Let us first detail the geometry of the problem. The two-point correlation function is 
defined as a function of the two correlated positions ri and r 2 , but on its expression (3.21) we 
realize that it can also be expressed as C^{ri,r 2 , r). Indeed, it is only a function of the shape 
of the triangle defined by the two sources and the observer, and this is fully characterized 
by the length of its three sides. This is because the global rotational invariance of the 
correlation function has absorbed three out of the six degrees of freedom. Furthermore, one 
can equivalently describe such triangle by two sides and an angle, see figure 2. Dehning (f> 
as the angle between the two directions {cos4> = i^), the correlation function can also be 
expressed in the form C^{ri,r 2 , (p) = C^{ri,r 2 ,r{ri,r 2 , (p))- The third side length r, and the 
length of the bisector of sources directions d (see figure 2 ) are obtained as 

r^{ri,r2,(p) = |r 2 - rip = (n + r2f sm^{(p/ 2 ) + (r2 - nf cos^{(p/ 2 ) . (3.22) 

In the plane-parallel limit, that is for (/> <C 1, then sin^(0/2) « (p"^/A, and cos^{(p/2) « 1 —0^/8. 
Therefore, 

« (r2 - nf + ^(ri 6rir2 -b rl)(p^ . (3.23) 

O 

We also assume that the two sources are far away compared to their separation, that is if 
r <C d, then d « ri « r 2 . We then approximate /3i ~ /32 f5{d) and Gi ~ G 2 ^ G{d) that we 

simply note /3 and G. From eqs. (3.23) we get immediately dr/d(/> Ri (r^-|- 6 rir 2 -|-r 2 )<^/( 8 r) Ri 
d‘^<p/r, from which we deduce 



Instead of using the set of parameters (r, d, (p) to express the correlation function in the 
plane-parallel limit, it proves useful to keep r and d, but to use 9 (or /rur = cos 6) dehned as 
the angle between r and the bisector (see figure 2). In the plane-parallel limit, it is related 
to (p by dcp/r Ri sind and cos 6 Ri |r 2 — ri|/r. 



Figure 2 . Representation of the geometry of the problem, ri, r2 and r form a triangle, r can also 
be expressed in terms of ri, r2 and the angle (p. 6 is the angle determined by the bisector and r. 
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We are now ready to examine the correlation function given by eq. (3.21) in the plane- 
parallel limit. The term linear in the differential operator is expressed in terms of jo(A:r) 
and jQ{kr)/{kr). Using (A.15) this can be simplified as 


2 


joikr{(j))) = --k^ [jo{kr)PQ{nAr) + 32{kr)P2{ndr)] (3.25) 


where the Pi are Legendre polynomials whose explicit forms are given in eqs. (A.23). For the 
term quadratic we obtain contributions from j^\ j^\ Jq, and Jq, where j^\kr) means 
the n-th derivative of the spherical Bessel function. Using the relations (A. 15) and (A. 16), 
this contribution can be cast as 
2 


1 /I 


16 


—jo(fcr)Po(lfdr) - l^j 2 {kr)P 2 {lJ,dr) + l^j4{kr)P4{fidr) 


15 


21 


35-' 


^ jo{kr{^)) = k 

(3.26) 

Inserting eq. (3.25) and eq. (3.26) into eq. (3.21), we finally obtain that the plane-parallel limit 
of the two-point correlation function in configuration space depends only on the difference 
of positions between the sources r. More precisely, there is an axisymmetry around the 
common line of sight of the sources d (that can be taken as the z-axis), so it depends on 
r, the distance between the two sources, and on its orientation with respect to the common 
line of sight which is given by /r, but not on an azimuthal angle. Indeed, it takes the general 
form 

?pp(d,r) = ^^p(r) = C'^p(ri,r2) = ^ ^f\r)Pe{fidr), (3.27) 

^=0,2,4 

where the only non-vanishing coefficients ^^(r) of this expansion are directly obtained from 
eqs. (3.21), (3.25), and (3.26) as 


^^(0 = (l + H[](r), if\r) = - (^/3 + ^“(r), e:’{r) = , 


35 

(3.28) 


and with the convenient definition 


k'^Ak 

27r2 


[kr) '^ji{kr)V{k) 


(3.29) 


Eq. (3.27) is the Kaiser formula in configuration space [22]. 

3.3 Redshift-space distortions in Fourier space 
3.3.1 General expression 

In this section, we investigate the Fourier conjugate of the objects given in (3.13) and (3.17). 
We first note that for a function of k • r we can relate operators acting on the r dependence 
to operators acting on the k dependence. Indeed for scalar functions of k • r. 


§^/(k-r) = ^/(k-r), 


= ^/(k-r) 


L?/(k.r) = L?/(k.r) 


Using these results, we find an equivalent expression for eq. (3.13) which is 


(3.30) 


(f^(r) = (1 -f /3(r))(5(r) -h 


/3(r)G(r) J 


= {1 + I3{r))6{r) + p{r)G{r) 


d^k ^(k)L? 

(27r)^/2 j-2 

d^k 6ik) , 

(27r)3/2 (A:r)2 k 


Jkr 


,ikr 


(3.31) 
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Once this step has been taken, it is then straightforward to take the Fourier transform to find 
a relation between the Fourier components of 5^ and 6. In principle we must only use the 
property that the Fourier transform of a product is the convolution of the Fourier transform. 
However, our goal here is to understand the geometrical structures induced by RSD effects 
in Fonrier space, and it will be much more transparent to consider the unrealistic case where 
j3{r) and G{r) can be treated as constants. In that case we simply replace j3{r) by j3 and 
G{r) by 1. We obtain, that 


5^(k) = (l + /3)5(k)+/3 
Furthermore, recalling that 


d^p 

p2 


5(p)L? 


gi(p-k)r 

(27r)3 


gi(p-k)-i 

(27r)3 


1 


47r|k — p| ’ 

we can rewrite this result in a more compact form with an integral on a kernel as 


(3.32) 


(3.33) 


d^p 


5^(k) = (l + /3)<5(k)+/3 / :^5(p)/C(k,p), /C(k,p)=L2^^ 


I 


47r|k — p| 


z/fcp = k • p . 

(3.34) 

Correlating a pair of fields using eq. (3.34) and the statistical properties (3.14) of the 
(I(k), we get 


C"(ki,k2) = (l + /3)2p(fei)4(ki-k2) + /3(l + /3) 


+/3^y d^p^^/C(ki,p)/C(p,k2) 


^^/C(ki,k2) + ^^/C(k2,ki) 


kl 


(3.35) 


The kernel K, deserves a closer analysis. An explicit form can be obtained using the expansion 


1 


k 


Ik - k'l 




h 
I 


(3.36) 


where fc< = min{/c, k''\ and A:> = max{/i:, k''\. We then find 


/C(k,k') =L 




1 


dvrjk — k'l 


dvr 


E 


K 




ki 


(3.37) 

because of the property (A.21). The coefficients of the expansion of the kernel JC in spherical 
harmonics are then easily found to be 


ICi{k,k') ^ ^ /C(k,k') = Y,l^eik,k')Yem{mUk '). (3.38) 

2^+1 tr 

The most important identities satisfied by JCi are presented in appendix B. This same kernel 
was already introduced by [12]. 

The two point correlation of the density fields in Fourier space is also the Fourier trans¬ 
form of the two point correlation function in confignration space. This means that eq. (3.35) 
should be related to eq. (3.21) by a Fourier transformation as 

C^(ri,r2) = ^ I d3kid3k2C^(ki,k2)e'*^i-^i-‘*^^-^= . (3.39) 
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In order to check this explicitly, we must proceed as follows. First, from (3.35) we use 
(3.37) to decompose the angular dependence in spherical harmonics. After taking the inverse 
Fourier transform, we need to also decompose the exponentials in spherical waves using the 
expansion (A.8). The angular integrals can be performed easily using the orthonormality 
of spherical harmonics. The final angular dependence can be simplified using the addition 
theorem for Legendre polynomials (A.24). As for the remaining radial integral, it can be 
simplified using the relation (B.8). The result is finally recast in the form (3.21) if we use the 
addition property (A.12) is employed to contract the spherical Bessel functions into jo(A:r). 


3.3.2 Alternative expression of the correlation function 

It is convenient to define another kernel to recast eq. (3.35) in a simpler form. Given that 
— 47 r|r-r'| Green’s function for the Laplace equation in three dimensions, i.e.. 




(3.40) 


and the decomposition (3.12) of the Laplacian, we define a kernel with the radial operator 
in Fourier space as 




1 


^47r|ki - ]^2 

With this kernel, the correlation in Fourier space reads simply 


= /C(ki,k2) + A;25^(ki-k2) 


(3.41) 


C"(ki,k2) = P(fei)<53,(ki-k2)+/3 


:^AA(ki,k2) + i^AA(k2,ki) 


kl 


+/3^y d^P^^^AA(ki,p)AA(p,k2). 


(3.42) 


This could have been obtained directly if we had started with eq. (3.11) and not from 
eq. (3.13) as chosen in the previous section. 


3.3.3 Plane-parallel limit in Fourier space 

As emphasized in section 2.2, the Fourier mode associated with the difference of positions is 
the median Fourier mode, and the Fourier mode associated with the median position is the 
difference of the Fourier modes, as can be seen on eq. (2.3). This is paramount to understand 
the structure of the Fourier space. It is thus natural to use the variables 


k = ^ , Ak = kl - k2 , (3.43) 

to parametrize the correlation functions in Fourier space. Defining = k • Ak and z^i 2 = 
kl • k 2 , they are related through 

kf = + {Ak/2)‘^ + k (A/c)/rkA (3.44a) 

kl = k‘^ + {Ak/2f - k (AA:)/xkA (3.44b) 

kik 2 h'i 2 = k^ — {Ak/2)‘^ (3.44c) 

^ 1 ^ 2(1 - ^ 12 ) = fc^(AA;)^(l - //kA) • (3.44d) 
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Since the plane-parallel limit in configuration space is the limit where r d, then in 
Fourier space it corresponds to |Ak| 3> k. Following this logic, we need to find the behaviour 
of eq. (3.42) under this limit. We find in this section that it is of the form 

Cp%(Ak,k) ^ C'^p(ki,k2) = dn{Ak)cP{k) + —^ Cf\k)Pe{f ,,,^). (3.45) 

' ' £=2,4 


The kernel AA(k, k') appearing in eq. (3.42) can be computed explicitly. However this 
requires to deal with the singular point k = k'. We thus need to use the relation (see eq. 
(3.46) of [23]) 


Ok 


where the superscript PV refers to the principal value. Using the geometrical relations (3.44) 
yields 


47rAA(ki,k2) 


2|k|2 


P2{lJ^kA) + 


1 

2|Ak| 


+ 47r^^(5D(Ak). 


(3.47) 


The first two terms are the regular part, and the last term is the distributional component of 
the kernel. This latter term contributes necessarily to the lowest order, that is to the plane- 
parallel limit, since it is necessarily diagonal in Fourier space. By dimensional analysis, the 
plane-parallel limit of the first two terms is necessarily proportional to l/|Ak|^, and this is 
another way to realize that the expansion must take the form (3.45). The plane-parallel limit 
of the kernel AAk(Ak) is then found to be 


47rAA(ki,k2) 


2k‘^ 

|A]^ 


^^(TfkA) + 


3 


<5D(Ak). 


(3.48) 


Using this limit, the plane-parallel limit of the Fourier space correlation (3.42) is 


C'p(Ak.k) = P(fc)«D(Ak) +2/!f^Vk(Ak) 

+ I i^Vk(Ap)Vk(Ak - Ap). (3.49) 


It is not expressed in the desired form (3.45) yet, due to the last integral. However if we 
notice that this is a convolution (defined with the (27r)^/^ factor) of a function with itself, 
then we can directly say that its Fourier transform is the product of the individual Fourier 
transforms. The inverse Fourier transform of the plane-parallel kernel is just 

+ Ia^=) = Afi . (3.50) 


where we used the Rayleigh formula to expand the exponential, and Weber integrals (see 
appendix A. 4) for the radial integration on spherical Bessel functions. Finally, the coefficients 
of the plane-parallel limit of the Fourier space correlation function, as defined in eq. (3.45), 
are 


Cr(A)) = (l + ^/3 + i/l2)p(A)), cf (A;)= (-4/3-y/l2^P(A)), C,fHk) = ^^'^V{k). 

(3.51) 
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This result is not the Kaiser formula in Fourier space as originally derived in [2]. In fact 
in the plane-parallel limit of the correlation function in Fourier space given by eq. (3.45) 
whose coefficients are given in eq. (3.51), there are still mode couplings. However, the Kaiser 
equation holds for independent Fourier modes, and is diagonal. It is thus clear that a further 
approximation needs to be performed to obtain an equation involving uncoupled modes. In 
fact, in order to find the plane-parallel result (3.51), one had to consider the configuration 
space part of the plane-parallel kernel in eq. (3.50), using the median distance d rather than 
the difference of Fourier modes Ak. This means that in order to find the plane-parallel expan¬ 
sion in Fourier space, we had to go through a mixed configuration/Fourier space temporarily. 
We shall find that the Kaiser formula arises naturally in this mixed space. And as mentioned 
in section 2.3, once the coefficients of the double expansion in this mixed space are known, 
the coefficients of the double expansion in the Fourier space can be deduced extremely 
easily, since they are related by simple numerical factors. The next section is devoted to 
the computation of the plane-parallel approximation in this mixed spaced, allowing then to 
explain how the Kaiser formula arises in the subsequent section. 


3.3.4 Plane-parallel limit in mixed configuration/Fourier space 

The plane-parallel limit of the correlation function in the mixed space is obtained by inverse 
transforming the Fourier space correlation (3.49). With the help of the Fourier transformed 
plane-parallel Kernel (3.50), we get 

~ r rl^Ak 

Cp%(d,k) = J ^^Cp%(Ak,k)e‘^*^-^ (3.52) 

= (2^)372 • (3-53) 


The angular dependence of this expression can be decomposed onto Legendre polynomials in 
the form 

5p(d,k)= ^ Pf (fc)PKl^kd), (3.54) 

£=0,2,4 

with coefficients 


K\k) 


2 1 n 

1 + x/3 + -/3" 

3 5 


yjk) 

(27r)3/2 ’ 


H^\k) = 


Vf\k) 


8.2 nk) 

35^ (27r)3/2 • 




nk) 

(27r)3/2 ’ 


(3.55) 


Alternatively, one could also obtain the plane-parallel limit in the mixed space by taking 
the Fourier transform of the configuration space correlation function (3.27) in that same limit, 
that is with 


gp(d,k) 


d^r 

(27r)3/2 


Cpp(d,r)e 


—ikr 


(3.56) 


Expanding the exponential in Legendre polynomials with the Rayleigh expansion (A.8), and 
then performing the integral on the polar angle using the orthogonality relations of Legendre 
polynomials (A.26), we get 


gp(d,k)= ^ 


£=0,2,4 


vf--/' 


drr^j£(A;r)^)'’^(r) 


( 0 ), 


Pi{.hy.d) 


(3.57) 
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The object inside squared brackets in this expression is the Hankel transform of the 
Inserting eqs. (3.28) in (3.57), and using the orthogonality of spherical Bessel functions (A.13), 
we can check that 

gp(d,k) = gp(d,k), (3.58) 

meaning that we can obtain correlation functions in this mixed configuration/Fourier space in 
two different, but equivalent ways. Either we Fourier transform ^^(d, r) on the second variable 
obtaining ^^(d, k), or we take and inverse Fourier transform on ((^(Ak, k) on the first variable 
obtaining C^(d, k). These two objects are intermediate steps between full configuration or 
full Fourier space quantities, and we have explicitly verified here their equivalence in the 
plane-parallel limit. 



Figure 3. Relations between ,^^(d,r), C'^(Ak, k), ^^(d,k), and C^(d,k). The hat denotes Fourier 
transform on the second variable, a tilde denotes (inverse) Fourier transform on the hrst variable. 
Arrows pointing to the right indicate Fourier transforms while arrows pointing to the left indicate 
inverse Fourier transforms. The objects inside the box are equivalent. 


3.3.5 Recovering the original Kaiser formula 

The first method to obtain the Kaiser formula consists in making the substitution /ikd 
^kz in tbe plane-parallel limit (3.53) of the correlation function expressed in the mixed 
configuration/Fourier space, i.e. we consider that the geometrical dependence is a function 
of the angle determined by k and the z axis. This particular z direction is chosen as being 
some average direction of the survey considered. With this extra approximation, we can take 
the inverse Fourier transform to obtain 

Cp%(Ak,k) « CLis«(Ak,k) = SD{Ak)V{k){l + (3pil,f = SoiAk) 

£=0,2,4 

(3.59) 

which is now the Kaiser limit as originally derived [2]. The main consequence of the Kaiser 
approximation, is that the correlation function is now perfectly diagonal, as it would be for 
a homogeneous distribution, but it still keeps a directional dependence. 

However, it is not a well defined object in Fourier space, since the angle /ikz makes 
reference to the direction of the z axis, defined in configuration space, and for which there is 
no unambiguous definition. The expression (3.53) now gives a clear meaning to eq. (3.59). 
Indeed, only the power-spectrum at a given distance d is well defined, and it has by definition a 
mixed dependence on variables in configuration and Fourier spaces. The crossed conjugation 
of variables in configuration and Fourier spaces allows to use geometrical arguments to define 
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a plane-parallel limit inside the framework of Fourier space correlations, but still referring to 
the position space. 

Since we have also shown in the previous section that the mixed space expression can be 
obtained from the configuration space correlation function, the second method to obtain the 
Kaiser formula consists in performing the replacement /ird A*rz in the expression (3.27) of 
^pp(d, r). Again the direction z is ambiguously defined, and is usually some average direction 
in the survey. In that case the dependence in d drops out, that is we use the approximation 


gp(d,r)»S,,..,(r)= E 

£=0,2,4 

When going to the full Fourier space we get 


(27r)3/2 J (27r)3/2 

Jd Ak f d^i 


J (27r)3/2 

= (27r)3/25^(Ak) 


^Kaiser 


(27r)3/2 


(r) 


(27r)3/2 


(3.60) 


(3.61) 


and with the expansion (3.60), we also recover the Kaiser limit (3.59). 


3.4 Multipole decompositions and angular correlations 

When calculating the effect of RSD on the correlation function in Fourier space, we derived 
that C^(ki, k 2 ) is given in terms of the kernel /C(ki, k 2 ) by eq. (3.35), and also that the kernel 
fC admits a decomposition in a basis of spherical harmonics as presented in eq. (3.38). It is 
natural to ask if C'^(ki,k 2 ) inherits the decomposition property of the kernel JC. We shall 
show that this is precisely the case, and after defining the multipoles in the next section, we 
derive in section 3.4.2 the coefficients Cf{k\, ^ 2 ) of the multipoles correlations. Furthermore, 
we will show that the kernel JC is also present on the expression of the distorted matter 
density field in configuration space, and we obtain the multipoles correlations C^{ri,r 2 ) in 
section 3.4.3. 


3.4.1 Angular multipoles and Hankel transformations 

We decompose the distorted and underlying matter density fields on spherical harmonics 
basis as 

,5(k) = Y,Sim{k)YeM , 5^(k) = ^<5|^(A:)Y,^(k), (3.62) 

£m im 

^em{r)Yim{r ), -5''(r) = ^ 5|^(r)Y£^(f). (3.63) 

im im 

From the statistic of the underlying density field (3.14), we find that the associated multipoles 
statistics is simply 

{Simirn'm'ik')) = . (3.64) 

Furthermore, using the Rayleigh expansion (A.8), we find that the multipoles in configuration 
and Fourier spaces are related by a Hankel transformation 

Sem{k) = J drr'^ji{kr)dimir), 6^rn{r) = J dkk‘^ji{kr)6im{k) (3.65) 


- 16 - 








with similar relations between the distorted fields multipoles 5l^{k) and 


3.4.2 Correlations of multipoles on Fourier space 

From eq. (3.34) and eq. (3.37) we find that the Fourier space multipoles are related by 

5l^{k) = (! + /?) 6em{k) + pj dk'JCeik, k')6em{k') . (3.66) 

Because of the global rotational invariance, the correlation of multipoles is necessarily diag¬ 
onal, even with the RSD effects. In particular, the correlations of Fourier space multipoles 
are of the form 

)) — ^U'^mm'Cl{k,k') . (3.67) 

The C^{k, k') are related to the Fourier space correlation function thanks to 

C"(ki,k2) = ^ J](2£ + l)C|(fci,fe2)P£(ki -ks) = J^C|(fci,fe2)Y£^(ki)Y;^(k2). (3.68) 

i im 


Inserting eq. (3.37) into eq. (3.35) we obtain, after integrating angular dependences, 


C^(ki, k2) = (1 + /3)' Soiki - k2)^%^ + /3 (1 + /3) 


kl 


V{ki) ^ V{k2) 


kl 


kl J 


V{p) 


-h/3^ / dp— j^iCe{ki,p)JCi{p,k2) 


p. 


^^^5{ki - k2) + + 


kl 


kl 


kl J 


+13'^ j dp^^J\fi{ki,p)Ni{p,k2), 


where we defined a new kernel 

Me{ki,k2) = lCi{ki,k2) + (^z)(fci, /C 2 ) • 


(3.69) 


(3.70) 


Note also that eq. (3.69) could be obtained directly from (3.66) using the statistics (3.64). 


3.4.3 Correlations of configuration space multipoles 

We have obtained the relation (3.13) from physical arguments, and derived its Fourier coun¬ 
terpart (3.34). We can, however, use eqs. (3.13), (3.20) and (3.40) to write 

6^{r) = {1 + I3{r))6{r) + J d^r'6{r')JC{r,r'). (3.71) 

It is striking that the same kernel appears, both in the configuration and Fourier spaces 
relations between the distorted and underlying matter density fields. From eq. (3.71) we find 
that the configuration space multipoles are related thanks to 

= (1 + / 3 {r)) Sem{r) + J dr'r'‘^}Ci{r, r')Sem{r') ■ (3.72) 

Alternatively, in the case where /3(r) is approximated by a constant function /3 as we did 
when analyzing the Fourier space structure, this relation could have been obtained using 
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a Hankel transformation (3.65) on eq. (3.66). This would only require to use eq. (B.7) to 
express }Ci{r, r'). 

Thanks to global rotational invariance, the correlation of configuration space multipoles 
is necessarily diagonal and thus takes the form 


The Cf{r,r') are related to the configuration space correlation by 

C^(ri,r2) = ^C|(ri,r2)Y£^(fi)Y^^(r2). 


(3.73) 


(3.74) 


im 


Using (A. 12) in (3.21) the angular dependence can be easily integrated, and this leads to 


C|(ri,r2) = GiG2 J k^dk^V{k) 


l-^O 

^2 


1 _ 

k? 


je{kri)je{kr2) 


(3.75) 


This agrees with the C^s given in [15] once restricted to RSD effects in wide angles. 

Using (A.7), it is easy to see that (3.75) can also be written under the alternative form 


C’|(ri,r2) = GiG2 J k^dk^P{k) 


1 + /3i - Mkn) [1 + /32 - /32-^J Mkr2) • 

(3.76) 

Note that in the case where /3i and f32 can be approximated by constants (and Gi = G 2 = 1), 
then G|(ri, r 2 ) could also be obtained from G|(/ci, /C 2 ) given by eq. (3.69) thanks to a double 
Hankel transformation, if we use eqs. (B.7) and (B.8). Alternatively, it is possible to obtain 
C^{ri,r 2 ) directly from We need only to take the spherical harmonics decomposition 

of eq. (3.16) to obtain 


J dkk‘^dim{k) 


1-^0 


je{kr). 


(3.77) 


Equation (3.75) is then obtained easily if we use the statistics (3.64). 


3.5 Representation of all derivation paths 

We can perform a diagrammatic representation of all relations obtained so far as in figure 4, 
so as to summarize the results obtained. 


• In section 3.1 we have obtained the general expression (3.13) of the distorted field 5^(r) 
from the underlying field due to RSD effects. 

• This allowed to compute the correlation function in configuration space G^(ri, r 2 ) (one 
southbound step in the diagram) in eqs. (3.17) and (3.21). 

• In section 3.2.2, the plane-parallel limit C'pp(ri,r 2 ) = Cpp(r) has then been obtained 
(another southbound step in the diagram), in eq. (3.27). 

• In section 3.3, we obtained the general expression for the density field in Fourier space 
(i^(k) in eq. (3.34) (one step in the north-east direction from (5^(r)). 

• The correlation function in Fourier space was then derived in eq. (3.35) (one southbound 
step). 
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Figure 4. Diagram representing the relations among quantities in configuration space, Fourier space, 
and their counterparts in harmonic space. We represent the various ways in which these objects are 
related, indicating invertible relations by double headed arrows. Plane-parallel limits are indicated 
with simple arrows. 


• The Fourier space plane-parallel limit was derived in section 3.3.3 (bottom line of the 
diagram) but it was found in Secs. 3.3.4 and 3.3.5 that we need to rely on a mixed 
configuration/Fourier space to be able to define correctly a Kaiser limit. 


• The angular multipoles were obtained in configuration and Fourier spaces in eqs. (3.72) 
and (3.66) respectively (one step in the north-west direction from (5^(r) and (i^(k)). 


• and Sf^{k) are related by Hankel transformations (3.65). 


• The multipoles correlations were obtained from the correlation functions in config¬ 
uration and Fourier spaces in eqs. (3.69) and (3.76) (closing the cube in the diagram 
from the bottom side), and they are related by Hankel transformations. 


• Furthermore, we explained that they can also be obtained by correlating and 

(5|^(A:), thus closing the cube from the north-west side. 
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4 Wide angle effects and plane-parallel corrections 


The remainder of this article is dedicated to the analysis of the corrections that should be 
added to the plane-parallel correlators to recover their general expressions, so as to take 
into account wide angle effects. We will thus explore in more details the bottom part of 
the diagram in figure 4. After expressing the general correlation function in a much more 
convenient form in the next section, we obtain the plane-parallel corrections in configuration 
space in section 4.2. We then deduce the corrections in the mixed configuration/Fourier 
space in section 4.3, and then we explain that this is enough to obtain the plane-parallel 
corrections in the full Fourier space in section 4.4. 


4.1 Improved correlation function in configuration space 

The expression (3.21) is general as it includes all the wide angle effects in configuration space. 
We now explicit the differential operators and o (noting 1^12 as v in this section for 
simplicity) which appear in /[jo(A:r)] so as to obtain an expression for which the corrections 


to the plane-parallel limit are more easily obtained. Dehning 

= (-rir 2 )" , hn = f jn{kr), (4.1) 

\ r dr / V ^ / 

we can show that the operators take the form 

= (1 - iy2)p2 _ , (4.2) 

O = (1 - l/2)2p4 _ _ (-4 3) 

The differential operators P” are closely related to those appearing in the expressions (A. 11 ) 
for the ji. Therefore, 

^Ikikr) = (1 - i'‘^)h2 - 2uhi (4.4) 

o L^jo(/cr) = (1 — — 8 ; 2(1 — z/^)/i 3 -|- (14z^^ — 6 )/i 2 -|- ii'hi. (4.5) 


Using these results, we find that 


mkr)] = (1 + A)(l +ft)/.o + [(1 - - 2^h,] 

+ ^kr!)^kr2p ^ ^ 

Using eqs. (A.14) it is then immediate to obtain an expression involving only jo, j 2 , and jh 
I[jo{kr)] = 


hikr) 

+j2{kr 


(1 + /3i)(1 + /32)- 


212 / ri 


2i 2 ( ri 


r2 


/32(1 + /3i) + -/3i(1 + /32) +1302 

r2 ri 


U0 - 6' 
15 


r2 


—4“ —/32(1 + /3i) 4-/3i(l + /32) ] + {1 — 12 ) 


r2 


ri 


2 n f ri(32{l + /3i) -t- r|/3i(l -k f32) 


o o I ^ /I 2\''’U2 280 — 12 

+0h ( ^ )-^ H 


+ji{kr) (302 


(1 - 0 ) 


^2^2 
2\2'1'2 


8 




rir2 1A:0 — 6 


rpA y rp‘2 


35 


o a 4i^r^ r jo(fcr) j 2 {kr) 

^ ^3rir2 (A;r)2 (/cr)^ 


(4.7) 
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We prove in appendix C that (4.7) coincides with the two point correlation function given in 
[1]. We have already obtained in section 3.2.2 that in the plane-parallel limit the geometrical 
terms multiplying Jq, j 2 and jh reduce respectively to terms proportional to the Legendre 
polynomials Pq, P 2 , and P 4 , leading to the Kaiser formula in configuration space. We will now 
examine the corrections to this plane-parallel limit and show that each of these geometrical 
terms is a linear combination of Legendre polynomials. 

4.2 Wide-angle effects corrections in configuration space 
4.2.1 Geometry of correlation functions 

The geometry of the RSD effects depends on the shape of the triangle whose sides are ri, 
r 2 and r. However there are alternative parametrizations for this triangle. A possible choice 
of parameterization which is examined in appendix C is to use two angular variables for the 
sources directions, and the distance r between them. 

However, ultimately we want to express our result as a power series in a small parameter 
which characterizes the departure from the plane-parallel limit. A first choice is to use the 
median parametrization 

{r,d,/rdr} with ^dr = d-r, ^ + ^2 

More generally, a family of choices can be parametrized by a parameter 0 < ?; < 1. We define 
d as the distance from the observer’s vertex to the point that divide the opposite side of the 
triangle, that is the line which connects the two sources, in the proportions v/{l — v). This 
family of choices is illustrated in figure 5. For each possible choice v, an angle 9^ is defined 
as the angle between r and d. In principle one should use the notation dy to specify that the 
average distance depends on the choice of the parametrization, but since no ambiguity can 
arise we always note it d. If ?; = 1/2 we recover the median parametrization. 



Figure 5. Given the triangle with sides ri, r 2 , r, we define a family of parametrizations r, d, 0^,, 
0 < V < 1. For each value of v, d is the length of the segment joining the vertex associated to the 
observer to the opposite side of the triangle, and 9y is the angle defined by the intersection of the 
lines d and r. 

Using that family of parametrizations, all the geometrical terms in eq. (4.7) can be 
expressed as a power series in r/d. In the plane-parallel approximation, that is for distant 
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sources, r and this parameter is indeed small. We further expand the angular dependence 
of each order in this expansion using a multipole decomposition on Legendre polynomials. 
The correlation function is thus expanded in a double sum as 

OO OO OO 

e(d, r) ^ C*(d - (1 - o)r, d + »r) = ( 5 )h"'”'(d. r) = E O” 

n=0 n=0 i=0 

(4.9) 

In the following we consider three parametrizations. The simplest choice is asymmetric and 
corresponds to u = 0. In that case the distance to one of the sources is used as an average 
distance since d = r 2 . We will show that this choice is not optimal. Then we consider the 
median choice u = 1/2 for which d = |ri + r2|/2. It is the choice described in section 2 and 
it arises naturally given the structure of the Fourier transform. Finally we also consider the 
choice in which d is the length of the line that bisects the angle (p as it might be the easiest to 
implement observationally. For simplicity, we ignore in the remainder the variations of the 
bias /S(ri) and /I(r 2 ) and we assume a constant bias (3. Similarly we neglect the variations of 
the transfer fonction G{r) and set G = 1. This is will simplify considerably the comparison 
with the results formulated in Fourier space and help emphasizing the geometrical aspects of 
the RSD effects. For each case we expand the operator I[jQ[kr)] in powers of r/d, restricting 
to second order corrections, as 

^bo(fcr)] = I^^^Mkr)] + ^I^^''[jo{kr)] + (^) ■ (4.10) 

Once inserted in eq. (4.7) it allows to obtain the coefficients Ci^\r) of our correlation function 
expansion (4.9). 


4.2.2 Asymmetric parametrization (v = 0) 

From basic triangle geometry, we obtain 



(4.11) 


where we abbreviate fJ-dr = k-- We first expand these relations in powers of r/d = r We 
then need to insert them in the geometrical coefficients of eq. (4.7). The relevant terms are 
collected in Table 1. 

The zeroth order correction, corresponding to the plane-parallel limit is 

/(°)[jo(A:r)] = ^1 -h jo(^r)Po(li) - {^13 + 32{kr)P2{k) + ^ 13 "^ ji{kr)Pi{k) ■ 

(4.12) 

Once inserted in eq. (4.7) we recover the plane-parallel result of eqs. (3.28). The non-vanishing 
coefficients in this plane-parallel limit are those of eqs. (3.28). 

The correction linear r/d is 




(Pi(f‘) - P3(f‘))+k§(AM - P,A) 


hikr) 


,2 16 


+/3 ^ [-P^ik) + Pbik)] k{kr) 

63 


(4.13) 
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Table 1. Expression of the geometrical coefficients in (4.7) when d and 6 are defined in the asymmetric 
case r; = 0. 


Following the same method, the non-vanishing coefficients (r) are 

= (i/3 + |/3") =§{>■). eiPir) = - (i/3 + |/3=) =§('■) - . 


= ^/3"=JM, 

where the H™(r) are integrals on the power spectrum already dehned in eq. (3.29). 
The quadratic corrections are 


(4.14) 




P ( F^o(/^) - lP2{^l) ]-0^[ ^Poif^) + 0-P2{f^) 


9 


9 


15 


joikr) 


+ 


/3(4p„(,, + Hp,(,)_|p,(,) 


( — Po{^J.) - —P2if^) + —P00) 

^ \ 35 147 245 


j2{kr) 


+ ^Pa(.) ) 74(.) 


+/3\ 


Mkr) ^ j2{kr) 


3 [ {krY {krY^ 


Poik^) 


(4.15) 
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Table 2. Expression of the geometrical coefficients in (4.7) when d and 6 are defined according to 
the median parametrization v = 1/2. 


(2) (2) 

Following again the same method, one should be able to read the Q . At this order, only Q , 
(2) (2) (2) 

would be non-vanishing. We do not report them since in this parametrization 
the second order corrections are subdominant. 

We remark that in this asymmetric parametrization the first order correction brings 
odd multipoles into play, as already observed by [24, 25]. For the second order correction 
there are only even multipoles, but beyond i = 0, i = 2, and £ = A corrections, there is a 
i = 6 term. In general higher order corrections bring higher multipoles to the scene. 


4.2.3 Median parametrization (v = 1/2) 

After some algebra, we find for this parametrization 


rf = (f + (r/2)^ + dr/Xdr (4.16a) 

r2=(f + ir/2)‘^ — drfifir (4.16b) 

rir 2 1^12 = dP - pr/2)^ (4.16c) 

rfrf(l - “ d-lr) ■ (4.16d) 


From eqs. (4.16) expanded up to second order, we obtain the terms needed which are gathered 
in Table 2. The lowest order results are the same as it does not depend on the choice of 
parametrization. The first order corrections vanish in that median parametrization since 
vanishes identically. We thus report only the second order which is 
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/3 n^0(//) - 7T^2(/i) - /3M -Po(//) + T7^2(/x) 


15 


jo{kr) 


+ 




11 


23 


48 


-/3 + r—P2{^^) + 7777 ^ 4 ( 1 ^) 


35 


147 


245 


j2{kr) 


16 


+/3 -7717^o(/i) - 7777^2(/i) + 7777T7^4(/i) + 7777^6(1^) 74(1^) 


+/3' 


105 


hikr 


{krY {krY 


245 ' 2695' 


231 


(4.17) 


The coefficients can be read from eq. (4.17), using eq. (A.14) to express j 2 {kr)/{krY in 
terms of jo, j 2 and jh, and the non-vanishing ones read 

= (?/3 - (I 7 ) Hg(r) + (7/3 - ^/S^) =«(r) + t/3^=g(.), (4,18) 

= TE^^) "St--) - (§'’ + f 

(?\r) = - (i/3 + ^/3^) Hg(r) + 2i5'9'=S(--). 




(4.20) 

(4.21) 


Because of the simplicity of the relations (4.16), we can explicitly verify that all the 
terms on the left column of Table 2 are expanded only with even powers oi r/d and p, to 
all orders. It implies obviously that = 0 when n is odd, and in particular that the 
first order corrections vanish. Furthermore, since only even powers of p appear in the series 
expansions at all orders, the parity of Legendre multipoles implies that only even is appear 
in the decomposition (4.9). Finally, we verify from the structure of the terms on Table 2 with 
the rules (4.16) that, for a given even order n in the series expansion, the highest possible 
power of the angular dependence is It follows that for a given n, the sum in i in 

eq. (4.9) does not extend to infinity, but is limited by imax = n + 4. 


4.2.4 Angle bisector 

If we take d to be the length of the line that bisects the angle </>, v will depend on the sides 
of the triangle. With the help of two basic triangle geometry results known as Stewart’s and 
angle bisector theorems, we relate the variables 


r2 



d? + vr"^ + r? 



\ — V V 

n r2 ’ 


(4.22) 


and then determine the geometric factors of eq. (4.7) in Table 3. As for the median parametriza- 
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Table 3. Expression of the geometrical coefficients in (4.7) when d and 6 are defined according to 
the bisector parametrization. 


tion, the first order corrections vanish identically. The second order corrections are then 
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7i^o(h) - - r 3To(/i) + 


+ 
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15 


29 


13 


29 


24 


-r + Tj=P2{fJ^) + 7777^4(h) 


35 


147 


245 


j2{kr) 


16 


( ~T(^Poik') + 7 ;;Tr-P 2 (/t) + 7777-Rt(h) J 4 (/^) 


105 


735 


245' 


jo{kr) 




jo{kr) j 2 {kr) 


Poif^). 


3 [ (kr)^ {krY 

It follows from (4.23) that the non-vanishing coefficients Q {r) are expressed as 


(4.23) 


if\r) 








4 

Is' 


- { -/3 + -/32 ) H0(r) - 




V 63 


147 


' 2 (^) 245 ^ 


2^0(,) 




(4.24) 

(4.25) 

(4.26) 


We observe that in this parametrization only corrections to the monopole, the quadrupole, 
and the hexadecapole are introduced. The fact that different multipoles appear in different 
parametrizations is expected as it is related to the natural sensitivity to the choice of the 
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origin in multipole expansions. In fact, we can verify also in this case that only corrections 
for n even appear, and that only even multipoles contribute but, for given n, imax = n + 2 
in the expansion (4.9). 

To summarize, in configuration space wide angle effects are encoded in corrections linear 
in r/d if the parametrization used is not symmetric enough, e.g. for the case u = 0. The 
largest effect in that case is the dipole correction {i = 1) which could be used to constrain /3. 
However for a balanced parametrization such as the median or the bisector parametrizations, 
the corrections are quadratic in r/d. In all cases, the only integrals on the power spectrum 
that need to be performed up to second order are 

/ dtp f dk 

-i^V{k)j(,{kr) with £ = 0,2,4, and El{r) = J ^^^V{k)jo{kr). 

(4.27) 


4.3 Wide-angle corrections in mixed confignration/Fourier space 

In the plane-parallel limit, we have already emphasized in section 3.3.4 that the standard 
results appear in the mixed configuration/Fourier space and not in the Fourier space. In that 
case we need to Fourier transform the dependence in r while keeping the dependence in the 
average distance d. This leads to define a power spectrum at a distance d. The dependence 
on d in the configuration space correlation function (4.9) appears through its norm d and 
through the angle /idr = d-r. We then find that in the mixed space, the natural perturbative 
expansion is in powers of l/{kd) as 


r(d,k) 


r] T* 

(27r)3/2 ^ 


(d,r) 




^Akd) 

n=0 


with an angular multipoles expansion given by 


(4.28) 


£=0 

with /Xkd = k • d. Equation (4.29) is obtained by expanding the exponential in spherical 
waves and using the addition theorem for Legendre polynomials. Given the global rotational 
invariance, we find that from the apparent six degrees of freedom of ^^(d,k), only three 
remain as it depends only on d, k and ^kd- 

At lowest order, that is considering the term n = 0 in (4.28), we get 


I Adr{kr)'^ji{kr)^^"'\r) 






nk) 

(27r)3/2 

Vjk) 

(27r)3/2 


1 + + 7/5^ -Po(hkd) + 

3 5 


+ UA P 2 M + 

3 7 35 


(1 -|- (dn 


kdt 


(4.30) 


thus recovering the plane-parallel result (3.54). 

The wide angle corrections can also be computed for each type of parametrization. 


Asymmetric parametrization (v = d) 














The first order correction is 




J p2dpP(p)|(^i/3 + ^/32 )X2i-i(p,A:)Pi(/ikd) 


+ 


^2 3-2(P,fc) + ^/32X4 3_3 (p,A:) 


-P3(/^kd) 


16 

+ ^/^^2i54_4(p, A:)P5(/Zkd) 
bd 


(4.31) 


where the integrals Xpg„ are defined in appendix F.l and listed in F.2. After integration 
by parts we obtain 




(27r)3/2 


V{k) ( ^/32P5(/ikd) + ) ^3(Mkd) 


+ ^ + §/3') ^i(/^kd)^ + kV'{k) - ^/32p5(^j^^) 

+ “ 45^^) PsifJ-kd) + -Pi(l^kd) ) • (4.32) 


We observe that this contribution to the power spectrum is purely imaginary, as ob¬ 
served by [26], and is expressed in terms of V{k) and V'{k). Furthermore, the coef- 
hcients V^^\k) can then be read directly from it. As a check, we perform the same 
analysis for u = 0 using cylindrical coordinates in appendix D. 

At second order, we obtain 


?^^^^(A;,/ikd) = J p'^dpVip) 


^/3- )Xooi(p,A:) 


+ ( l20o{p,k) - ^13^1,0_,{p,k) 
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+ Q W (^^oooCp, k) +Io2-iip,k)) 
3 p^ 


PoiPkd) + 
12 


g/3 + y^/3^ ) Po2o{p,k) 


^/3+ ^/3'J P22-l{p,k) - ^/32X42_2(p,fc) 
[2,2 ^ . 96 ^ 776 2 


-P2(Mkd) 


-^/5^2:46_4(p, A;)P6(l^kd) I • 


/3 X44-3(p, A:) 


Pi{p. 


kdj 


(4.33) 


Similarly the integrals Xpqn are defined in appendix F.l and listed in F.2. We do not 
report the final result at second order since it is subdominant in that case. 

• Median parametrization (u = l/2j 
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The first order correction vanishes, and at second order we get 




^/3- ) Xooi(p,A:) 


+ ( ^/3 - 2 : 200 (P, A:) - A:) 


4/52 

+ 0 ^ (2:ooo(p, k) +Xo 2 -i(p, A:)) 
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-Po(Pkd) 
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~*~245^ 2 : 42-2 (p, A:) 


-P 2 (Pkd) + 
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A(pkd) - 2:46_4 (p, A:)P6(Pkd) ? .(4.34) 


Following the same method, after integrations by parts we obtain the second order 
correction 
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. (4.35) 


/Q'f 

The coefficients y {k) are read directly from this result. We remark that it is expressed 
in terms of V{k), V'{k), and V"{k). In general at n-th order it will involve the ra¬ 
th derivative of V{k). As a check, we perform the same analysis for v = 1/2 using 
cylindrical coordinates in appendix D. 


Bisector parametrization 



















The first order correction also vanishes, and at second order we get 
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(4.36) 


Following the same method we are able to Hnd that the second order corrections are 
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(4.37) 


i‘2') 

from which the corresponding {k) are easily read. 

As we can observe from section F.2, the integrals Xpqn appearing in (4.33), (4.34), and 
(4.36) are all expressed in terms of Dirac delta function and its first and second derivatives, 
except for Z 4 o_i(p, A:) and Io 2 -i{p, k), which also involve Heaviside step functions. The 
numerical coefficients, however, are such that the contributions from the Heaviside step 
functions always cancel out. 


4.4 Discussion and Fonrier space correlations 

We notice that among the coefficients of the plane-parallel expansion, some are always 

vanishing. Indeed they only take non-vanishing values if (. and n are either both odd or both 
even. For instance at lowest order the non-vanishing coefficients are ^o^^(r), and 

'^ 4 ^^(^)- At first order in the asymmetric case, the non-vanishing coefficients are ^(^^(r), 
and '^ 5 ^^(r) whereas at second order we have only C^\r), ^ 4 ^^(r) and 

In fact, we have shown at the end of section 4.2.3 that if we use the median parametriza- 
tion the only non-vanishing coefficients are those for which I and n are both even at all orders. 
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Including second-order corrections would lead to an expression which is correct up to fourth- 
order corrections. The number of terms to be included is limited since I < n -|- 4. The 
coefficients in the mixed space follow obviously the same structure since they are related by 
the Hankel-type transformation (4.29). The inverse transformation takes the form 





k^dk{kr) ’‘ji[kT)'pf^\k) 


(4.38) 


Given that the coefficients in the mixed space are expressed only in terms of the matter power 
spectrum and its derivatives, then provided it has nice convergence properties at /cd —)• 0 and 
fed —>• oo, the transformations between configuration space and mixed space coefficients are 
in general not worrisome. 

However, if we want to build an expansion of the correlation function in the full Fourier 
space, it is not obvious at first sight that it is well defined. Indeed, following the logic that 
the Fourier mode associated with the difference of points is the average Fourier mode, and 
the Fourier modes associated with the average of points is the difference of Fourier modes, 
we are lead to consider that the small parameter in the plane-parallel expansion in Fourier 
space is |Ak|/fe. We are thus tempted to look for an expansion of the type 

C(Ak,k) = —^ ^ with /xkA^k-Ak, (4.39) 

the prefactor l/|Akp being introduced from dimensional arguments. Performing a Fourier 
transform on ^(d, k), and using the Rayleigh formula (A.8) to expand the exponential, would 
lead to Hankel-type relations between the coefficients in mixed spaced and Fourier space in 
the form 


C!['\k) = 'P^''\k)]J^i-iY j (4.40) 

Vf'\k) = Ce'^\k))J^i^ J x'^~^ji{x)dx. (4.41) 


From the Weber integrals of appendix A.4, we already see that this would not be defined 
if i and n were not both odd or both even, that is if they would not have the same parity, 
given that one would possibly encounter the poles of the F functions, e.g. for n > £ 3. 

Since we have shown that £ and n always have the same parity, such problem never arises. A 
second problem remains for £ = n = 0 though, since it involves F(0). If we inspect the plane- 
parallel limit obtained in eqs. (3.45) and (3.51), we notice however that there is a monopole 
multiplying a Dirac delta function, and it differs from the attempted expansion (4.39). Hence, 
we should consider instead an expansion of the form (2.5) where the factor dvr has only been 
introduced for convenience. Using eqs. (A.18), we then recover that the relation between the 
Fourier space and mixed space multipoles is given by eqs. (2.8). It can be checked that at 
lowest order, that is in the plane-parallel limit, the Fourier space coefficients (3.51) and the 
mixed space coefficients (3.55) are indeed related by eqs. (2.8). Since we have already shown 
in section 3.3.5 that the Kaiser formula naturally arises in the mixed space, we conclude 
that one should always rely on a mixed space expansion rather than on a full Fourier space 
expansion when apprehending the RSD effects in Fourier space. 


- 31 






5 Conclusion 


The wide angle effects on RSD have a rich geometrical structure, both in configuration space 
and Fourier space. This richness is highly simplified in the plane-parallel limit where we 
get an expression for the two-point function that does not depend on the distance from the 
observer to the pair of galaxies and, on the other side of the coin, independence of Fourier 
modes in the Kaiser formula. 

Our initial goal was to detail how such tremendous simplifications occur, and to provide 
intermediate formulas between the full wide angle treatment and the plane-parallel limits. 
Instead of using the formulation introduced by [1] where a functional basis built from the 
monopole of the tensor product of three sets of angular functions is used, we derived a very 
simple expression for the general correlation function, given in eq. (3.21), in terms of an 
integral on the product of the power spectrum V{k) with the action of differential operators 
on the spherical Bessel function jo{kr), where r encodes all the information about the triangle 
formed by the observer and the pair of galaxies. Beyond its simplicity, this formulation of 
the problem proved to be very powerful to demonstrate relations between configuration and 
Fourier spaces. 

The Fourier conjugate of eq. (3.21), which is eq. (3.35), paves the way for unveiling 
the geometrical structure of the approximations leading the Kaiser formula. Indeed, the 
presence of the kernel JC both in configuration and Fourier spaces umbilically connects both 
formulations, and its (singular) structure leads to the expected plane-parallel limit. However, 
we argued that a mixed configuration/Fourier space has to be employed, since it is only in 
this space that the plane-parallel limit is meaningful. Indeed, if one looks at the Kaiser 
formula, the dependence on the angle between a Fourier mode and the line of sight cannot 
be well defined neither in Fourier space nor in configuration space. It is the structure of the 
plane-parallel limit in the mixed configuration/Fourier space that gives a precise meaning to 
this phantasmal dependence. 

Going back to the general expression for the two-point function given in eq. (3.21), we 
worked out the action of the differential operators on the spherical Bessel function to ob¬ 
tain an explicit formula in terms of lengths and angles given by eq. (4.7). We first observed 
that this expression coincides with the general result presented by [1] (as demonstrated on 
appendix C), and then proposed a perturbative representation to this general result, as pre¬ 
sented in eq. (4.9). On this perturbative solution we looked at contributions classified by 
powers of r/d, where r is the distance between two galaxies on a pair and d is a mean distance 
from the observer to the pair. In the plane-parallel limit this ratio goes to zero, and only 
the term n = 0 contributes on the sum. The coefficients are the plane-parallel limit 
ones, and are only non-vanishing for .^ = 0,2,4. Departing from the plane-parallel limit, we 
found the r/d corrections, and the corresponding corrections in l/{kd) in the mixed config¬ 
uration/Fourier space. However, these corrections depend on the parametrization used to 
express the geometry of the triangle in eq. (4.7). If we choose an asymmetric parametriza¬ 
tion we can show that dipolar terms appear at first order, as shown in eq. (4.14). Such a 
parametrization could be used if one is interested in probing the values of (3. 

As for symmetric parametrizations, we have two natural options for the distance d, 
namely the median and the length of the angle bisector. We show that if these parametriza¬ 
tions are used, the first correction to the plane-parallel limit is quadratic in r/d, and therefore 
the full wide angle corrections are reduced even for large coverage surveys, as already found 
numerically by [6, 7, 27]. If one includes second order corrections, the parametrization of the 
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angle bisector proves to be the most immediate extension to the plane-parallel limit, since 
only £ = 0,2,4 contribute in eq. (4.9). 

Finally, we have shown that it is possible to define a plane-parallel expansion in the 
full Fourier space. However, its coefficients are related to the mixed space coefficients by 
simple numerical factors and one should instead consider the mixed space expansion as the 
spectrum at a given distance is the natural quantity to extend the notion of power spectrum. 

Acknowledgments 

The authors thank Lado Samushia, Marc Manera and Raul Abramo for discussions about 
redshift-space distortions. We are indebted towards G. Faye for his help about the structure 
of the RSD kernels. This work has been done within the Labex ILP (reference ANR-10- 
LABX-63) part of the Idex SUPER, and received financial state aid managed by the Agence 
Nationale de la Recherche, as part of the programme Investissements d’avenir under the 
reference ANR-ll-IDEX-0004-02, and also under the reference ANR-12-BS05-0002. P. R. 
also thanks FAPESP for financial support. 


A Special functions 

We collect the some properties of Bessel functions and Legendre polynomials in this appendix, 
since they are frequently employed throughout the main text. 

A.l Bessel functions of the first kind 

Bessel functions of the first kind are solutions of Bessel’s differential equation 


iA 

X dx 




^2 


Ji{x) 


(A.l) 


regular at r = 0. For a negative index, the dehnition is extended through J_£ = (—l)”j£. 
They can be obtained from the integral representation 


= dBe 

J2tt 


irsine-Ue _ t ■\l 


1 


= {-iY-l d9e 

J27r 


ir cos Q—\iQ 


(A.2) 


This is particularly useful to write the Fourier series of Indeed from this integral 

representation, we obtain 


cos Q 






£=—oo 

Two recurrence relations satished by Jf are of particular interest 


j;(x) = 


Ji-i(x) - Ji+i(x) Ji{x) j£+i(x)j£_i(x) 


2 x 2i 

Finally, for all £ € N, the Bessel functions can be generated from Jq from 


(A.3) 


(A.4) 


Jii{x) = x^ 



(A.5) 
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A.2 Spherical Bessel functions 

The spherical Bessel functions are related to Bessel functions of the first kind by 


k{x) = 



and they are the solutions regular at the origin of the differential equation 


J_A 

J.2 




V r2 


k{r). 


(A.6) 


(A.7) 


The spherical Bessel functions appear naturally in a spherical geometry when a plane wave 
is expanded into spherical harmonics or Legendre polynomials. Indeed, this is given by the 
Rayleigh expansion which is 


h"-*" = 4vr ^ \^jkkr)Y^m{x)Y}^{k) = '^{2i + l)i^je{kr)Pe{k • f) (A.8) 

im I 

where the are the spherical harmonics and the Pi are the Legendre polynomials. Two 
recurrence relations satished by ji{x) are particularly useful 

jdx) ^ k+i{x) +ji-i{x) 

X ~ 2^+l 

{2i + l)j'i{x) = iji-i{x) - {£ + l)ji+i{x). 

Spherical Bessel functions for any n € N can be obtained from jQ{x) thanks to 

k{x) = {-l^x^ , joix) = . 

yxdx/ X 

If r = r 2 —ri, then e^^'^ = e‘*^ R 2 -ri) Rayleigh expansion (A.8) and the addition theorem 
for Legendre polynomials (A.24) imply that 


(A.9) 

(A.IO) 

(All) 


jo{kr) = '^{2£ + l)ji{kri)ji{kr 2 )Pi{ri ■ f 2 ) (A.12) 

£ 


where r, ri, and r 2 must form a triangle (r^ = rf + — 2r 1 X 2 X 1 ■ f 2 ). Note that this is a 

special case of Gegenbauer addition theorem (eq. (3) on section 11.4 of [28]). 

The spherical Bessel functions satisfy the orthogonality relation 

j drr‘^ji{ar)ji{br) = ^<^ 0(0 - b). (A.13) 

A.3 Useful relations among spherical Bessel functions 


The relations involving spherical Bessel functions of low order and their derivatives which we 
have used in this article are 

jl{x) j 2 {x)+jo{x) j 2 {x) l./N, 2.,, 23(.x) j4{x)+j2{x) 


15- 


21 -^ 


35“ 


7 


= -^ihix) + j2ix)) , foix) = -^jo(®) + ^j2ix) 


(A.14) 

(A.15) 


( 0 ) 

= \jo{x) + ^j 2 {x) - ^k{x), jo\x) = ^jo(a:) - ^j 2 {x) + ^k{x) ■ (A. 16) 
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A.4 Weber integrals 

From eq. (1) on section 13.24 of [28] we obtain the following Weber integrals 

r[(p + £ + i)/ 2 ] 

Jo ^ ’ r(^/ 2 -p /2 + 1 ) 


(A.17) 


In particular, if I and n are both odd or both even, such that we can define s by £ = 2+n + 2s, 
we obtain 


- / ji{x)x "dx 


I je{x)x'^ Mx = 

The main Weber integrals used in this article are 

J j 2 ix)x~^dx = i , 

j j2(x)x^dx = 3, 


-1 


2”+"(n + s)! 
(2s+ 3)!! 


dx = 


2 

TT 


2 f 2 
- / J 2 (x)x dx = 


15 

15 

Y 


(A.18) 

(A.19) 

(A.20) 


A.5 Legendre polynomials 

Legendre polynomials are solutions regular at the origin of 


dfi 






where — 1 < /i < 1. They obey the parity relation 

The more relevant Legendre polynomials in this work are 

^ 3/i2_l 35Y-30/i2 + 3 

Poip,) = 1 , P2iii) = ---, P^ilJ.) = -r- 


(A.21) 


(A.22) 


(A.23) 


Legendre polynomials satisfy an addition theorem, expressed in terms of spherical harmonics 

^r((f 1 ■h) = Yl i)YL(f2). (a.24) 

m 

They are also used to dehne a closure relation for the angular part of a Dirac delta function 
with 

''^^ ^^^Y(ri-f2) = ^D(ri-r2). (A.25) 

For a function / of a variable p € [—1,1], using the orthogonality relation 


I' ^Pe{n)Pe{n)dn = 


(A.26) 


we can write a Legendre series decomposition as 


27 + 1 

= —~ / • (A.27) 

£=0 


- 35 - 















B The kernel K,i 


In order to obtain important identities satisfied by the kernel coefficients JC^, we first consider 
the integral 

— J drr‘^j'^{kr)f^{k'r)kk', (B.l) 

where the prime indicates differentiation with respect to the argument of the functions. We 
will compute this expression using two different methods. First we will use integration by 
parts, and then we will use recurrence relations. By means of integration by parts, we can 
show with the help of (A.7) that 


2 

vr 


drr‘^j'g^{kr)j[{k'r)kk' = doik — k') — —+ 1) 

vr 


drji{kr)ji{k'r ), 


(B.2) 


where we also have used (A.13). For the second method, we use eqs. (A.9), (A.10), and 
(A. 13) to obtain 

j drr‘^j[{kr)j[{k’r)kk' = Snik - k') - j <^rrji{kr)ji_i{k'r) 

+k j drrje{k'r)ji_i{kr)^ . (B.3) 

Each of the integrals on the right hand side of (B.3) belongs to a class of discontinuous 
integrals of Bessel functions. The equation (8) on section 13.42 of [28] can be used to show 
that 




= 


A (AV 

K k' \k') 


k' < k 
ii k' = k 
\i k' > k 


(B.4) 


Hence, we find that 


^ j drr‘^j[{kr)j[{k’r)kk' = doik - k') - + i) 

where /c< = min{/i:,A:'} and /c> = max{A:,A:'}, and using the dehnition (3.38) it is simply 
expressed as 

— J drr'^j[{kr)j[{k'r)kk' = doik — k') + lCf,{k, k'). (B.6) 

Comparing (B.2) and (B.6), we conclude that 


lCi{k,k') 


TT 


drjf,{kr)ji{k'r). 


(B.7) 


We finally obtain a useful relation derived from the orthogonality relation (B.7) 

[ dkk‘^}Ce{k,k')je{kr) = -^^—^^je{k'r ). 


(B.8) 
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C Comparison with [1] 

We shall now demonstrate that the two point correlation function given in (4.7) coincides 
with the expression given in [ 1 ] when evaluated in the case of constant bias (/3i = /32 = /3). 
The two-point correlation function in configuration space is written in that reference as 


= ^l dk^Pik) 

+bni n 2 sin(ni 4>i) sin(n 2 4>2) 


flni 712 cos(ni (pi) cos(n 2 ( 1 ) 2 ) 

. 711,112=0,1,2 


(C.l) 


where the non-vanishing coefficients are given by 


aoo = 3o{kr) - Q/3 + j2{kr) + k{kr) 

ao2 = 020 = + ^ld‘^k{kr) 

022 = i^k{kr) - ^j2{kr) -h /3^ 


^22 = 


^jo{kr) - ^j2{kr) - ^k{kr^ 
(fcr) 5^ {kr) 


Ooi —- 

92 


2/3 + -/3' 

5 


^o 2 \k{kr) I 2 h{kr) 


{kr) 


5^ {kr) 


(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C. 6 ) 

(C.7) 


Oil — 


012 = 


1 f 4 jo(fcr) _ 8 j 2 {kr) 

9192 V3 {kr)"^ 3 (A:r)2 

1 f2ji{kr) 8jkkr) 


/3' 


\5 {kr) 5 {kr) 

5^2 = — f + lk{kr) 


gi \5 {kr) 5 {kr) 






hii =- 


‘^jo{kr) ^ 4: j2{kr) | ^3 


9192 V3 {kr)"^ 3 (fcr)^ 

1 (2ji{kr) 3j3{kr)\ 3 

02 I —-1 - X -- /, X ) P 


92 \5 {kr) 5 {kr) 


92 \5 {kr) 5 {kr) 


(C.8) 

(C.9) 

(C.IO) 


with 


91 


sm{(p2) 


sin((/72 - Pi) 


92 


sm{4>i) 


sin((^2 


- <Pi) 


(C.ll) 


Our strategy is to sum explicitly all the terms in eq. (C.l) using some properties of 
spherical Bessel functions and trigonometric identities to obtain eq. (4.7). For this sake, we 
consider the partial sums 


Si= ^ anin 2 Cos(ni(/)i)cos(n 2 (/> 2 )+ 6 nin 2 sin(ni(/)i)sin(n 2 (/> 2 ), i = 0,1,2 (C.12) 

()ll,7l2)GAi 
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with Ao = {(0,0)}, Ai = {(0,1), (1,0), (1,2), (2,1), (1,1)1, and Aa = {(0,2), (2,0), (2,2)}. 
With the help of eqs. (A.9), (A. 14) and the trigonometric relations 


cos (f)l + cos (/>l COs(2(^2) cos 4>2 + COs(2(^i) cos (/>2 


gi 

cos (/>1 cos (f >2 
gi 52 Vsin (/)2 ' sin(/)i 

cos(2(;i)i) cos (/>2 cos cos( 2 (/) 2 ) 


= -4 


52 


sm4>i sm(j)2\ , , , ^ 

= 2-1 ^^^ ) COS ((/>2 - (pi) 


— —2 + 2sin^((/)2 — (f>i) -\- 


52 51 

sin(/>i sin( 2 (/> 2 ) sin(2(/>i) sin (/>2 


51 


52 


= -2sm‘^{(p2 - 4>i) 


sm^{4>2 - 4>i) 

sin (pi sin (p 2 


sin (pi sin p 2 
sin p 2 sin pi 


cos pi cos p 2 , 

(C.13) 

(C.14) 

COs{p2 - Pi) , 


(C.15) 

(C.16) 


we can express Si as 


^4 sin ^((/>2 - pi)cos{p 2 - pi) f 5 o(fcr) j 2 (fcr) \ ^ 

^ 3 sin (;i)i sin (/)2 \ (fer)^ (fcr)^ ) 

2 /si^ cos ((/>2 - (/)i)(jo(A:r) + j2(A;r))/3(l + /3) 

3 \smp 2 sin01 J 


+ -^{jo{kr) + j2{kr))P{l + /3) + sin2(02 - Pi) 


4 6 2 

-Tjo{kr) - -j2{kr) - —jiikr) 

5 7 35 


/3' 


The trigonometric identities 

cos(20i) cos( 202) = —1 + (cos(20i) + cos(202)) + 4sin^ 01 sin^ 02 , (C.17) 

sin(20i) sin(202) = 2 — 2sin^(02 — 0i) — (cos(20i) +cos(202) — 4sin^ 0i sin^ 02(,C.18) 
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allow us to write S 2 as 


15 


-52 = - - — sin ^((/)2 - 0i) 13 jo{kr) + 1 - -(cos(2(/)i) + cos( 202 )) I3j2{kr) 


+ “TTT - T7(cos(2(/)i) + cos( 2(/>2)) + — sin^(02 - 4>i) P'^hikr) 


21 14 

12 8 


21 


35 35 


7 


+ ( sin^((/)2 - 4 >i) + ^(cos(2(/)i) + cos( 2(/)2)) + sin"* sin^ 412 ) ^^jAikr). 


Finally, using the identities in Table 4, we can check that (/>i, 42 ) = ^ / (Ikk"^P{k)[So-\- 
5i + 52 ] is precisely eq. (4.7). 


D Cylindric method 


We first write the configuration space correlation function as in [1]. For simplicity we neglect 
the variations of /3(ri) and f3{'f2) and consider a constant bias /3. We also neglect the variations 
of the transfer function and set G = 1). We get 


r(d,r)^C^(ri,r2) = 

r d^k 


d^k 

(27r)3 


1-^0 
^2 '^^1 


(2vr) 


1 _ (k- ri)^ 




k^ 


} - —O 

^ J .2 *^^2 « 

2ik • £-2 (k • f2''^'' 

1 + ^ 


k?r2 


A;2 


ikr 


The operator which applies on the exponential is thus [21] 


0^1 ^2 — 


^ 2ik-fi ^ (k-fi)^ 


k'^ri 


fc2 


j ^ 2ik • h ^ (k-f2)2 


k4r2 




(D.l) 


(D.2) 


Like in the spherical method, we expand this operator in orders of plane-parallel perturbations 
as 

Or„ = 0<% + + (I)" OW,. (D.3) 

However, the difference is that we perform the integral on k using cylindrical coordinates 
instead of spherical coordinates. This implies that we need to decompose everything using 
these cylindrical coordinates. The natural choice for the azimuthal axis is the direction 
defined by d. The separation vector r = r 2 — ri is then decomposed into an axial part which 
is r^r with = d • r, and a polar part r_L whose length is just rcor where uj^ = 1 — that 
is 

r = /j,rd -|- r_L , r± = rcor , ry = rfir ■ (D.4) 

A general Fourier mode k is decomposed similarly in cylindric coordinates with an axial 
part which is along d, and a polar part which is orthogonal to this axial direction. The polar 
part has a length kujk with = 1 — but it is not necessarily aligned with r_L. That is 
we use 


k = fe/ifcd -|- k_i_, kj_ = kujk , ^11 = k^k > cosOk = r_i_ • k_i_ (D.5) 

k • r = /i;||r|| -|- k±r± cos 9k = kr{^rk-k + cos Ok) ■ (D.6) 

In order to perform the integration in cylindrical coordinates, we need to express all products 
of the type k • ri or k • r 2 in the operator (D.2) in terms of the radial and polar parts of k, 
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fi and r 2 - This depends on the choice of parametrization, that is on the choice of d which 
is given by the parameter v (see section 4.2.1). 

If we define as the angle between ri and d, and similarly for ^ 2 i then from basic 
trigonometry we get 


cos <I>i = d • f 1 = Y 1 — 
(1 — v)r 


r(l — v) 


ri 


UOr 


sin<I>i = 


-UJr = 


(1 — v){r/d)ujr 


ri 


i/l + (r/(i)2(l — v)^ — 2(r/d)(l — v)fj, 


cos <I '2 = d-f 2 = Wl — 


rv 

Sm <$2 = - OJr = 


Ur 


rv 
^2_ 
{r/d)vur 


^2 ^ + (r/(i)2ti2 + 2{r/d)vfi 


(D.7) 

(D.8) 

(D.9) 

(D.IO) 


Note that in the limit where the separation is much larger than the average distance d, the 
directions ni and n 2 are nearly d since <hi ~ <I >2 ~ 0, meaning that at lowest order k-fi = fij. 
and k • f 2 = fJ^k- Expanding in powers of r /d we find up to first order 

- ^ T - T 

k • ri =/Tfc - (1 - u)-a;fca;r cos^fc , k • f 2 =+ u-Wfc cos 0^ . (D.ll) 

If we consider the particular case of the median parametrization v = 1/2, we get up to second 
order 


k • ri = Hk- cos 4 + ^ {-d-k^l - 2/rr, Uk Ur cos Ok) (D.12) 


r 

J.2 

M2 


k • f 2 = /ifc + — Wfc Wr cos 9k + ^ OJk COS 9k) 


(D.13) 


Using these expansions, we can find the operator (D.2) order by order, in terms of the axial 
and polar parts. The lowest order is just the Kaiser operator 


— (1 + • 

At first order, if we choose the asymmetric parametrization (u = 0) we get 

= -2 Wfc/rfc(l + Pul) Ur cos 9k ■ 


(D.14) 


(D.15) 


In the median case {v = 1/2), there are no first order corrections, but the second order ones 
read 


^ [l^r^J-k{'i- + PA) - - P^J^l)uJr cos 9k] 


(D.16) 


+ ^ [-ulnl{l + Pill) - 2Ukllk{^ + Plll)llrUrCOs9k +^1(1 “ /3/r|) (cj^ COS 6»fe)^] . 

The correlation function is found from these operators by integration over 9k using the 
property (A.2) of Bessel functions. At lowest order, that is in the plane-parallel limit, we 
obtain 


rW(d,r) = 


(1 + Pill)" Jo(fc^r^)e''=ilbl . (D.17) 
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The corrections are defined using the expansion (4.9). At first order in the asymmetric case 
u = 0, we get 

= j + /3^fc)9fc^Jo(fc±r_L)h^ll’'ll . (D.18) 

At second order in the median parametrization (v = 1/2) there are no first order corrections 
and the second order corrections read 

f dk\\k±dkj_ /3f4o2 r, 2 <, 

+ UJk^lk{l + l3fil)dk,dk^ 

+ /3/ri)ri| Jo(A:^r^)e'"ril . (D.19) 

We can then check that the correlation functions found in eqs. (D.17) (D.18) and (D.19) 
match the ones obtained using the spherical method if we use the relation 

1'^ ^ r Jm PPif^k) = i^h{kr)pp{^ir) ■ (D.20) 

This general relation is simply found by expanding e'*^ ’’ either in spherical coordinates with 
eq. (A.8) or in cylindrical coordinates with eq. (A.3) and by integrating both expressions. 
Indeed, the integrand k±dk± in eq. (D.19) can be changed to kdk, and then the integration 
on /c|| can be performed using this formula, eventually leading to the results obtained with 
the spherical method in section 4.2. 

The spectrum at a given distance d is defined as in the spherical case (4.28), but the 
integral of the Fourier transform is done using cylindrical coordinates, leading to 

f'^(d,k) = J 27rdr||r_Ldr_LJo(p±r-_L)?^(d,r)e“'^ril , (D.21) 

where the integration over the polar angle in spherical coordinates has already been per¬ 
formed. We use the integrals of appendix E.2 in order to perform the integral on r±, and for 
the integral on ry we use the property 

J dry = 2^4">(fcy -py). (D.22) 

The hnal step consists in taking integrations by parts. 

With this method, we recover that the spectrum at lowest order is the Kaiser limit in 
Fourier space. For this we need only J^oo{p±,k±_) as we get after integrating over ry and r_L 

?<“>(«:,W) = + (D.23) 

= + (D.24) 

We can compute the spectrum at first order for the asymmetric case u = 0. For this we need 
only 77io(p±, k_i) defined in appendix E, and applying the same method as for the background 
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order, and after integration by parts, we get exactly the same expression as in eq. (4.32). 
Finally, in the median parametrization (n = 1/2), there are no first order corrections and the 
second order corrections are exactly the same as those of eq. (4.35). 

When deriving this result, the intermediate steps involve Jq{p±, k±), Jio(p±, k±), J 2 q{p±-, k±) 
and j7oi(px, kx_) of appendix E which are handled with integration by parts. 


E Integrals over products of Bessel functions 
E.l General Method 

We follow the method developed in the appendix of [29] to compute general integrals involving 
a product of Bessel functions. We start from the simplest orthogonality relation 

J xdxJo{ax)Jo{bx) = ^ ^ = jQ{a — h). (E.l) 

Then we define the following integrals 

Jnp = J xdxJn{ax)Jo{bx)x"'^‘^^ (E.2) 

xdx [d” Jo(ax)] Jo(6x)x”^^^ . (E.3) 

From the recurrence relations of the Bessel functions, we find that these integrals are also 
related by recurrence relations, which are 



i7n+l p(Q) ^) — 

1 n+1 (o, 6) — ^)] 

piP'i b) — da'J^npifl': b) ■ 

E.2 Integrals used in this article 

In this paper we have used only a few integrals of this type which are 
Jio{a,b) = daJoia^b) = ^ ^ 

J2o{a,b) = dadaJoia^b) = ^ ^ 

jQi{a,b) = -dadaJo{a,b) - ^daJo{a,b) + -^Jo{a,b) 

^ - b) _ S'pj{a-b) <5/,(a-6) 

b ab ba^ 


(E.4) 

(E.5) 

(E.6) 


(E.7) 

(E.8) 

(E.9) 


F Integrals over products of spherical Bessel functions 
F.l General method 

We follow a similar method for integrals of products of spherical Bessel functions. However, 
the class of functions that are required for this article is slightly more general. Let us dehne 
the integrals 

2 r°° 

Ipqn{a,b) = - drr2+2^+P+%(ar)jg(6r), (F.l) 

TT Jo 
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2 r°° 


drr'^j(,+p{ar)j^{br) ^ p2:£(a, 6) = X(£+p) ^ („/2_i_£_p/2) (a, ^) • (F.2) 

Obviously, we have Xpqn{a, b) = Zqpn{b, a). We will derive recursive relations using the recur¬ 
rence relations for the spherical Bessel functions. From 


i + l I 

= j£-i(x)-= -n+^ix) + -ji{x) (F.3) 

we get the following recursions 

2p_i q n+i{a, b) = -^^da [a^’+%gn(a, b)] (F.4) 

q—1 n+l(^5 ■db[W+%qn{a,b)] (F.5) 

Xp+l q n(a, b) = --^da [a~PIpqn{a, b)] (F.6) 

Ip q+i nia, b) = --^db [6"%q„(a, b)] (F.7) 

which can be generalized to give 

( 1 \ ^ 

[a~^Ipqn{a,b)\, (F.8) 


or in particular 


[a-n{a,b)] . 


(F.9) 


In order to use these recursions, we need to start from a known integral. The simplest starting 
point is just eq. (A. 13), which in the notation of this section is written as 


I{cL, 6) — Qlp(^(i, 6) — Xpp_p(a, 6) — 


Soia - b) 


(F.IO) 


Other starting points are necessary in this paper as we cannot reach all integrals from the 
previous starting point and the recursion relations. We find them by using eq. (8) of section 
13.42 of [28]. It leads to 


^-1 

ilp{a,b) (o, 6) qp+1 ~ 


(F.ll) 


where H is the Heaviside step functions which satisfies 77(0) = 1/2, H{x) = 1 if x > 0 and 
H{x) = 0 if X < 0. Then, the recurrence relation (A.9) implies that 


llp-i{a, b) = Zp+^p.^ _p{a, b) = {2p + l)^H(a - b) - 
Comparing (B.4), (B.6), and (B.7), we also conclude that 


(F.12) 


o2/(a, 7) — X££_(£_,_i)(a, 6) — 


1 


27 + 1 


Oj 

-^H{a-b) + -^H{b-a) 


(F.13) 
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F.2 Integrals used in this article 


With this method and starting from eq. (F.IO), we obtain for instance 

144-3(0,6) = (F'lS) 

2ooi(a,b) = — \daa^daT{a,b). (F.16) 

If we start instead from eq. (F.12) we also get 

l20-i{a,b) = ^ Ilia,b) (F.17) 

a-^X40-i(«,&) = a-^lX{a,b). (F.18) 


Using recursions similar to these two examples, we are able to compute the integrals needed 
in this article which read 


k‘^Iooo{k,p) = Snik-p) (F.19) 

fc%„4{fc,J.) = + (F,20) 

= sXiXX. _ /nik-p) _ j,) (p 21, 

k'‘l2a-i(k,p)^3X!iZll-SD(k-p) (F,22) 


»,=l40-l(t,rt = - 57^^ + _ s'Hk-p) (F.23) 


k^X22-i{k,p) 

k‘^Xo2oik,p) 

k^Xi2-2ik,p) 

k‘^X2A-2{k,p) 

k^Xi4^_2,{k,p) 

k^XiQ_^ik,p) 


= 8 


^ jMizi!) _ Mizrt _ (t _ p) 

Kk'P K p 

&Dik-p) , 6'jj{k-p) 


+ 5- 


+ 6Uk-p) 


15 

24 


p^ p 

boik-p) , J'jjik-p) _^f^{k-p) 


(kp) ^ k 
Soik-p) J^g^oik-P) 


P 


- ^oik-p) 


P 


+ 6'ljik — p) 


(F.24) 

(F.25) 

(F.26) 

(F.27) 

(F.28) 

(F.29) 
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k^Z2i-i{Kp) = -S'oik-p) (F.30) 

k‘^lA3-3{k,p) = _ s'^(^k-p) (F-31) 

k^l23-2{k,p) = + d'oik - p) (F.32) 

fc^X45-4(fc,p) = 6'j,{k - p). (F.33) 
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